Visualizing quantum mechianics in pliase space 
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We examine the visualization of quantum mechanics in phase space by means of the Wigner function and the Wigner 
function flow as a complementary approach to illustrating quantum mechanics in configuration space by wave func- 
tions. The Wigner function formalism resembles the mathematical language of classical mechanics of non-interacting 
particles. Thus, it allows a more direct comparison between classical and quantum dynamical features. 
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1. Introduction 

Quantum mechanics is a corner stone of modern physics and 
technology. Virtually all processes that take place on an 
atomar scale require a quantum mechanical description. For 
example, it is not possible to understand chemical reactions or 
the characteristics of solid states without a sound knowledge 
of quantum mechanics. Quantum mechanical effects are uti- 
lized in many technical devices ranging from transistors and 
Flash memory to tunneling microscopes and quantum cryp- 
tography, just to mention a few applications. 

Quantum effects, however, are not directly accessible to hu- 
man senses, the mathematical formulations of quantum me- 
chanics are abstract and its implications are often unintuitive 
in terms of classical physics. This poses a major challenge in 
teaching and learning the foundations of quantum mechanics 
and makes it difficult for students to form a mental picture of 
quantum dynamical processes. Therefore, teaching quantum 
mechanics may be supported by visualizing the mathemati- 
cal objects of quantum mechanics' mathematical language in 
order to sharpen the students' understanding of abstract con- 
cepts. In fact, many text books'^ communicate quantum me- 
chanics with the aid of visual means. 

Visualizing quantum dynamics is typically implemented by 
providing a pictorial representation of the quantum mechani- 
cal wave function in position space - or in momentum space''. 
In order to establish a relation between quantum dynamics and 
classical dynamics, the phase space seems to be more appro- 
priate because classical Hamiltonian mechanics is formulated 
in phase space. Thus, we will present in this contribution 
methods to visualize quantum dynamics in phase space. 

The remainder of this paper is organized as follows. In sec- 
tion [2] we briefly review how to depict quantum dynamics in 
position space or in momentum space using wave functions. 
In section [3] we give a short introduction to quantum mechan- 
ics in phase space and establish some connections between 
quantum mechanics and many-particle classical mechanics be- 
fore we visualize quantum dynamics in phase space in sec- 
tion E] 



2. Visualizing quantum dynamics in 
position space or momentum 
space 

A (one-dimensional) quantum mechanical position space 
wave function maps each point x on the position axis to a 
time dependent complex number *I'(x, t). Equivalently, one 
may consider the wave function 'i'ip, t) in momentum space, 
which is given by a Fourier transform of ^{x, t), viz. 

^^P^'^^^lk^^ J>I'(.,Oe--^/''d.. (1) 

Various methods to visualize complex-valued functions in po- 
sition space or momentum space have been devised. Figure [T] 
exemplifies three different ways of visualizing a wave func- 
tion. It shows a Gaussian wave packet 
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in position space with mass m, mean position x, mean momen- 
tum p, and width cr at time f = 0. 

Part a of Fig. [T] presents the complex wave function (O by 
plotting its real and its imaginary part as separate real-valued 
functions. Although a complex function is completely deter- 
mined by specifying its real and its imaginary part it is difficult 
to infer physically relevant information from these separated 
quantities. 

Visualizing the wave-function's squared modulus and its 
phase is more appropriate because the squared modulus of a 
wave function may be interpreted as a probability distribution 
and the change of its phase encodes information about the mo- 
mentum. Thus, we picture the wave function (|2]i in part b 
of Fig.[T]by small directional arrows. Each arrow's length is 
determined by the wave-function's modulus, its direction (an- 
gle between the arrow and the horizontal axis) by the wave- 
function's phase. The arrows in Fig.[T]b are sometimes called 
phasors. The method to portray a wave function by phasors 
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FIG. 1 : (Color online) Three different ways of visualizing the Gaus- 
sian wave packet (|2j at time f = with x = 1, p = 3, and cr = 1/2 
(dimensionless units with ft = 1 and m = 1); a) visualization by 
the wave-function's real part and imaginary part, b) visualization by 
the wave-function's modulus and its phase using arrows (phasors), c) 
visualization by plotting the wave-function's squared modulus and 
encoding the phase in a periodic color scheme. 

had been used, for example, by Feynman^i. It is beneficial in 
situations where advanced printing methods, color printing for 
example, are not available. However, it is difficult to picture 
fast oscillatory wave functions in this manner 

Another way of presenting a wave function by its modu- 
lus and its phase has been popularized by Thaller-ii^. This 
method is utilized in part c of Fig. [1] It shows the squared 
modulus of the wave function ^ plus its phase which is en- 
coded by the color between the graph of the squared modulus 
and the horizontal axis. Plotting the squared modulus allows 
for a direct read off of the probability distribution in position 
space. Note that the color scheme for the phase in part c of 



Fig. [U is 27r-periodic and contiguous. Encoding the phase 
as a color allows to plot even relatively strongly oscillating 
wave functions. It requires, however, high-quality color print- 
ing for good reproduction. Thaller'-- has also generalized this 
approach to multi-component wave functions. 

3. Quantum mechanics in pliase 
space 

3.1 . The Wigner function 

The Wigner distrubution function^^ was introduced by Eu- 
gene Paul Wigner to take into account quantum corrections 
to classical statistical mechanics in 1932. Today, it finds ap- 
pUcations in current research in various branches of physics 
ranging from quantum opticsiii and quantum chaos'- to nu- 
clear and solid state physics. The Wigner function allows 
for a mathematical description of (non-relativistic) quantum 
mechanics in phase space. But it is more than a sheer 
mathematical object, it may be reconstructed by a series of 
measurements^^— . The Wigner function formalism is math- 
ematically equivalent to the standard language of quantum 
mechanics in terms of wave functions and the Schrodinger 
equation but it is more closely related to classical mechanics 
because it formulates quantum mechanics in phase space as 
Hamiltonian mechanics does. 

The Wigner function of a (one-dimensional) pure state 
wave function *I'(x, t) is defined as 

\v{x, p,t)^-^ f '¥\x + s/2, 0'i'ix - s/2, t)e'''"'ds , (3) 
27rn, J 

where *P*(jc, t) denotes the complex conjugate of^'{x, t). Note 
that due to the symmetry of the integrand in (O, the Wigner 
function w(x, p, f) is always a real function. The expression 

w{x, p,t)^ f '^'•(p + s/2, tmp - s/2, t)e-''''^ds (4) 
2nn J 

is equivalent to (|3]l if ^(x, f) and "^{p, t) are linked via ([T]i. 
Further important mathematical properties^Ji of the Wigner 
function include the following relations. 

• If ^{x, t) is normalized to one then also w{x, p, t), viz. 

w{x,p,t)dpAx^ 1, (5) 
furthermore, the equation 

JJ w{x, p, tf dp = ^ (6) 
holds for all pure states '^{x, f). 
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Quantum mechanical probability densities in position 
space and momentum space may be obtained from the 
marginals of the Wigner function 

p(x, = mx, f)P = J w{x, p, t) dp , (7a) 

pip, t) = mp, t)f = J w{x, p, t) dx . (7b) 

The state overlap of two wave functions *i'i(x, f) and 
^2{x, t) in terms of their Wigner functions w\{x,p, t) and 
W2{x, p, t) reads 



%(X, tyi>2{x,t)dx : 



2nh 



p,t)w2(x,p,t)dpdx. (8) 



• The Wigner function obeys the reflection symmetries 

*P(x, f) 'i'*(x, t) => w(x, p, t) -> w{x, -p, t) , (9) 
T(x, t) '¥i-x, t) => wix, p, t) wi-x, -p, t). (10) 

• The Wigner function is Galilei-covariant, that is 
^{x,p,t) -^^'{x+y,p,t) 

w(x, p,f) ^ w(x + y, p,t) , (11) 

but not Lorentz-covariant. However, Lorentz-covariant 
generalizations of the Wigner function have been consid- 
ered in the literatursi^"— which may be appUed in relativis- 
tic quantum optics^l. 

• The Wigner function is not gauge invariant. The definition 
(O, however, may be generalized-- such that it becomes in- 
dependent of the gauge in the presence of a magnetic field. 

The Wigner function w(x, p, t) is called quasi probability 
function because it is real and normalized to one but not 
strictly non-negative. In fact, a sufficient and necessary con- 
dition that w{x, p, f) is strictly non-negative is that *F(x, t) is 
a minimum uncertainty wave packet, that is a Gaussian wave 
packet— . As a consequence of (|6j, the Wigner function can 
not be arbitrarily highly peaked, which is a reminiscence of 
the quantum mechanical uncertainty relation. 

The evolution equatio n^' "'^"^i^^ of the Wigner function may 
be obtained by calculating the Wigner-function's time deriva- 
tive. Taking into account the Schrodinger equation 

i/i— ^ = H>F(x, t) = + V{x, tmx, t) (12) 

at Lm ax^ 

for the wave function moving in the potential V{x, t) we get 



dw{x, p, f) 
~dt 



2^ J [{H^'ix, f))'P(x, ty - "Vix, t){H'V{x, t))*) e'l''" ds . 



After several algebraic manipulations that include expanding 
the potential V{x, t) into a Taylor series around x = and in- 
tegrating ( HjI ) by parts one finally finds the so-called quantum 
Liouville equation 



5w(x, p, t) 
Jt " 

2-J{2l+iy.\dp,,dxH dpHdx,, 



21+1 



H(x,p, t)w(x,p, f) , 



(14) 



which is a partial differential equation of infinite degree. The 
operator (37-^ - a^af") ^'^ * fT4b consists of dififerential op- 
erators that act on the Hamilton function H{x, p, t) (indicated 
by the index H) and differential operators that act on the 
Wigner function w{x, p, t) (indicated by the index w). The 
function H{x, p, t) in ( 114b is the classical Hamilton function 



H{x,p,t)^ ^ + Vix,t) 
2m 



(15) 



that corresponds to the Hamilton operator H of the 
Schrodinger equation il2\ . 

For systems where magnetic fields are absent as in ( 112b . the 
Hamilton function H(x, p, f) can be written as a sum of a ki- 
netic energy term K{p) - p^/{2m) that depends only on the 
momentum p and a potential energy term V(x, t) that depends 
only on the position x (and possibly the time t). In this case, 
the quantum Liouville equation (fT4l) may be written as a con- 
tinuity equation 



dw{x, p, t) 
Jt 



dx 
d_ 



■j(x,p,t), 



(16) 



where we have introduced the two-dimensional vector field 

^Mx,p,t) 



j(x,p,t) 



\ 1=0 



y (ih/2f (da a \-' av(x,t) . 
fi^aM)\[a^— - —.) -g]-w(x,p,t) 



(13) 



(17) 

called the Wigner function flow. 

Let Ox denote some operator function of the quantum me- 
chanical position operator and 0{x) the corresponding classi- 
cal function, then it follows from ( iTbb that 

{Ox){t) ^ {'¥{xj)\dx\'¥(x,t)) = J J wix,p,t)0{x)dpdx. 

(18) 

This means one can calculate the quantum mechanical expec- 
tation value (Ox) (0 of an operator Ox by averaging the clas- 
sical function 0{x) over the Wigner function w(x, p, t). For 
momentum dependent operators Op and their corresponding 
classical functions one finds using ( iTab a similar statement 

{0„) (f) = (^(x, t) I Op I ^^(x, f)) = JJ w(x, p, t)O(p) dp dx . 

(19) 
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For general operators Ox,p and classical functions 0{x, p) that 
depend on the position as well as on the momentum one can 
show^iSi that 

{0,,p)(t) ^ {^•{x,t)\6,,p\'V{xJ)) = 

w(x,p,t)0(x,p)dpdx (20) 



provided that the relation between the operator (9 v./j and the 
function 0{x, p) is established by 

d,,p = JJ (5(^, ny^^-'^^P^I^ A^An, (21) 

where 0{^, n) denotes the Fourier transform of 0{x, p) 

0(^,7r) = 0{x,p)e-'^'^*'"'''l''dxdp (22) 

and X the position operator and p the canonical momentum 
operator, respectively. 

The Wigner function we{x,p) of an energy eigen-state 
wave function of the Hamiltonian H{x,p) - p^/(2m) + V{x) 
is time independent and therefore it follows from ( I14l i that 

(p d V iih/2f d^'^'Vjx) 3^'+' ^ 

(23a) 

This equation, however, is not sufficient to specify the en- 
ergy eigen-value E, which is determined by the eigen-value 
equatiorJi 



£1_^^ y (i^/2)^' d^'V(x) d- ^ 
2m 8m5;c2+Zj {2l)\ 



we(x, p) 



dx'-' dp^'j 

Ewe{x,p). (23b) 

The two equations ( I23ab and ( I23bl l together are equivalent to 
the time independent Schrodinger eigen-value equation 



where H{x, p, t) denotes the Hamilton function that also gov- 
erns the motion of the ensemble's single particles. The equa- 
tion ( 126b may be written as a continuity equation ( fTST i by in- 
troducing the probabihty flow 



j(x,p,t) = 



£w(x,p,f) 
.SX^wix,p,t)) 



(27) 



2m dx- 



+ V{x, f)*P£(x, f) = E'¥e(x, t) . 



(24) 



The quantum Liouville equation (114b reduces in the limit 
?i ^ to the classical Liouville equation (l26T l. Note, how- 
ever, that also if the third spatial derivative of the potential 
V{x, f) and all its higher derivatives vanish the quantum Liou- 
ville equation (114b also reduces to the classical Liouville equa- 
tion (|26] |. This means that the quantum mechanical evolution 
of the Wigner function for a particle in a homogeneous con- 
stant potential or in a harmonic oscillator (or a combination 
of both) follows the same dynamical law as the probability 
density of an ensemble of classical particles in the same po- 
tential. Nevertheless, non-classical features may be present in 
the Wigner function. 

In summary, quantum theory in phase space in terms of the 
Wigner function features formal parallels to the classical the- 
ory for an ensemble of classical non-interacting particles. The 
state of both kinds of systems is completely characterized by 
a real-valued phase space distribution function w{x, p, t). The 
evolution of w{x, p, t) is determined by the classical Liouville 
equation (|26] | or the quantum Liouville equation (fT4l) . respec- 
tively. The quantum Liouville equation equals its classical 
counterpart plus some additional quantum terms. For classi- 
cal systems, w{x,p, t) may be any non-negative normalizable 
function. The quantum mechanical Wigner function, however, 
has to be compatible with the laws of quantum mechanics, 
e. g., to respect the uncertainty relation, and may be negative. 
Thus, proper quantum characteristics of a quantum dynami- 
cal process are governed the additional terms in the quantum 
Liouville equation or the negative parts of the Wigner func- 
tion. In fact, the negative parts of the Wigner function are a 
common measure for non-classicality'"*-'^'^'' -'^. 



3.2. Relations to classical mechanics 

The motion of an ensemble of classical non-interacting 
particles may be characterized by a probability distribu- 
tion w{x,p,t) in phase space.— The probability distribution 
w{x,p, f) determines how likely it is that a particle in the en- 
semble is at position x and has momentum p at time t. It 
allows to calculate expectation values of observables for the 
classical ensemble by an integral over the phase space 

<C,,p)(0 = JJ w{x,p,t)0{x,p)dpdx (25) 

likewise for a quantum system by (|20] i. The evolu- 
tion of w{x,p,f) is determined by the (classical) Liouville 

97 9S 

equatio n 1 

dw(x,p,t) dH(x, p,t) dw(x, p,f) dH(x, p,f) dw(x, p,f) 
dt dp dx dx dp 

(26) 



4. Visualizing quantum dynamics in 
pliase space 

After our short review of the Wigner-function's properties in 
the previous section, we are now prepared to visualize quan- 
tum mechanics in phase space. This may be accomplished 
by picturing the Wigner function itself or quantities that are 
derived from the Wigner function, e. g., the Wigner function 
flowdnll. 

In contrast to the wave function, the Wigner function is 
a real-valued quantity. It assigns just a single real value to 
each point in phase space, which simplifies the visualization 
of Wigner functions to some degree as compared to the visu- 
alization of wave functions, which assign two real values to 
each point in configuration space. The dimension of the phase 
space — which is twice the configuration-space's dimension — 
turns out to be a limiting factor when visualizing Wigner func- 
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FIG. 2: (Color online) Visualization of the Wigner function l l28t of a Gaussian wave packet by false color plots at time r = (left part) and 
t = \ (right part) with x = 1, p = 3, and cr = 1/2 (dimensionless units with h = \ and m =1). The panels above the false color plots show the 
probability distribution in position space l lTat while panels right next to the false color plots show the probability distribution in momentum 
space l lTbt . Values of theWi gner function are indicated by the color-bar at the left. 



tions. Therefore, we restrict ourselves in this contribution to 
systems with a two-dimensional phase space. Displaying two- 
dimensional real-valued functions is a standard task of scien- 
tific visualization and may be accomplished by various meth- 
ods, e.g., by two-dimensional false color plots (sometimes 
called heat maps), by contour line maps, which show lines of 
constant values of the Wigner function, or by (pseudo) three- 
dimensional surface plots. 

Systems with two or more spatial dimensions have a phase 
space with at least four dimensions which prohibits visualiz- 
ing the full Wigner function. In this case one may in order to 
reduce the dimensionality, for example, picture lower dimen- 
sional cuts of the Wigner function through the phase space 
or projections of the Wigner function onto lower dimensional 
surfaces. 



4.1 . A free wave packet 



H{x, p, f) - p^/{2m). Thus, the Wigner function of a free par- 
ticle satisfies the relation 



w(x, p, t) — w{x - pt/m, p, 0) . 



As a consequence of (|29] l, we find the equation 

p(x, ~ J" ^(^' P' t)'ip - J" — pt/m, p, 0) dp , 



(29) 



(30) 



which states that the probability distribution in position space 
p(x, t) at position x and time t is given by an integral over the 
Wigner function at time zero along phase space strips which 
are rotated against the p-axis by an angle of -arctanf/m. 
Thus, measuring the distribution p{x, t) of a function of 
the position x and the time t allows to reconstruct the full 
Wigner function w{x, p, 0) at time zero by means of an inverse 
Radon transformatio n "'-^ ' . This has also been demonstrated 
experimentall)^— . 



The Wigner function of the Gaussian wave packet Q is given 
by the strictly positive bivariate normal distribution 

1 / (x-pt/m-x)- 2crHp-pf 
w(.,p,r)=^exp^ — 

(28) 

which is illustrated in Fig. |2]by two false color plots for two 
different points in time, for f = in part a and for f = 1 in 
part b, respectively. At time f = 0, the Wigner function of 
the Gaussian wave packet is symmetric in phase space with 
symmetry axes parallel to the jc-axis and the p-axis. During 
the evolution of the wave packet, the Wigner function gets ro- 
tated and stretched along the x-axis in phase space. This cor- 
responds to the familiar wave packet broadening in position 
space, see also the probability distribution in position space 
dTal i as shown above the false color plots in Fig. |2] 

As we pointed out in section 13.21 the quantum Liouville 
equation ( 114b of a free quantum particle's Wigner function 
reduced to the classical Liouville equation i26i for an en- 
semble of free non-interaction particles with the Hamiltonian 



4.2. Characterizing eigen-states 

The Wigner function eigen-state of some time independent 
potential V{x) is determined by the two equations (12 3 al l and 
( |23bt . For the harmonic oscillator potential V{x) - maP'x^/l 
we find**'" the eigen-states 

w„{x,p) = 

^^'^^(-(^-^^'^^^^"(^(fe^^H' '"'^^ 



and eigen-energies 



E„ = hoj n 



(31b) 



with K - y/moj/h, n - 0, 1, . . . and L„{x) denoting the nth 
Laguerre polynomial. 

In Fig. [3] we depict the Wigner function of the harmonic 
oscillator eigen-state with energy E - ha){3 + 1 /2) together 
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with its Wigner function flow (fTTT i. Because the «th Laguerre 
polynomial has n positive roots, the Wigner function of the 
harmonic oscillator eigen-state oscillates and changes n times 
its sign. Note that despite the fact that the Wigner function 
( I31al ) is time-independent there is a permanent circular steady 
Wigner function flow in phase space as indicated by the arrows 
in Fig. [3] 

In section [3] we demonstrated that for harmonic potentials 
quantum dynamical effects are absent. This means, the dy- 
namics of the Wigner function follows the classical Liouville 
equation because all quantum terms in the quantum Liouville 
equation vanish for harmonic potentials. Also in equation ( |23] ) 
quantum terms vanish, which has consequences for the geom- 
etry of the Wigner function flow of the eigen-states of the har- 
monic oscillator Using equations ( fTTl ) and ( |23a| ) one can show 
that j{x, p) ■ Vw(x, p) - and, therefore, the Wigner func- 
tion flow is always tangential to the lines of constant Wigner 
function values. This is a feature that the harmonic oscillator 
eigen-states share with all time-independent classical phase 
space distributions. Also for classical distributions the flow 
( 127b is tangential to the lines of constant phase space distribu- 
tion. 

If the potential has non-quadratic terms then the quantum 
terms in ( I23l l will affect the structure of the Wigner eigen- 
functions and the Wigner function flow and the lines of con- 
stant Wigner function values are no longer tangential to each 
other Thus, the divergence between the direction of the 
Wigner function flow and the tangentials to the lines of con- 
stant Wigner function values may be interpreted as a measure 
how strong the non-classical terms in the quantum Liouville 
equation ( 114b affect the structure of the eigen-states. 

In Fig.|4]we show the Wigner function ground states of the 
anharmonic oscillator potential V{x) - marx^ /2 + ax'^ for var- 
ious values of the anharmonicity parameter a. For this figure 
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FIG. 3: (Color online) False color plot of the Wigner function of 
the harmonic oscillator eigen-function with energy E = 3 + 1/2 (di- 
mensionless units with h = 1, m = 1, and tj = 1). Gray solid lines 
indicate levels of constant Wigner function values while black an^ows 
indicate the steady Wigner function flow J17b . The arrows' length is 
proportional to the flow's magnitude. 



the Schrodinger ground state wave functions have been calcu- 
lated by imaginary time propagation----^^ and the Wigner func- 
tion was obtained by calculating (O numerically afterward. In 
Fig.|4]the anharmonicity parameter a grows from left to right 
starting with a = (part a) and increasing to a = 1/4 (part b) 
and a - 3/4 (part c). One clearly sees in Fig.|4]how the diver- 
gence between the direction of the Wigner function flow and 
the tangentials to the lines of constant Wigner function values 
grows with increasing a. 



4.3. Scattering by a potential barrier 

Scattering of a wave packed by a potential barrier is a well 
known textbook system to exemplify quantum effects. De- 
pending on the potential's shape and the initial condition of 
the wave packet some part of the wave packet may be reflected 
while the other part is transmitted; a process that has no anal- 
ogon in classical physics of point-like particles. 

The quantum dynamical scattering dynamics is commonly 
visualized by the wave function in position spaced and some- 
times in momentum space^. Figure |5] illustrates scattering of 
an initial Gaussian wave packet by a potential barrier in phase 
space at different points in time. The scattering potential bar- 
rier is located around x - and its shape is given by 

Vix) = yo(tanh(|)+ l)(tanh(-|) + l) , (32) 

see also Fig. |5]d. For this figure the Schrodinger wave func- 
tions was propagated by a Fourier split operator method-4 and 
the Wigner function was obtained by calculating ([3]) numeri- 
cally afterward. 

The system in Fig. |5] starts at time t - (Fig.|5]a) with an 
asymptotically free Gaussian wave packet. Here we consider 
a thin potential with height Vo - 64 and characteristic width 
A = 1/8. The wave packet's initial kinetic energy equals 
/ikin - 32. Accordingly to the laws of classical physics the 
kinetic energy is too small to cross the barrier. As the wave 
packet moves towards the potential barrier it starts to interact 
with the barrier. Due to the non-classical terms of the quan- 
tum Liouville equation the Wigner function starts to develop 
regions with negative values such that at about t - 1/2 sig- 
nificant interference patterns have emerged, see Fig. |5]b. At 
about f = 1 the wave packet has split into two parts that are 
traveling into opposite directions in position space, a reflected 
part and a transmitted part that has tunneled trough the po- 
tential barrier In phase space this corresponds to two isolated 
wave packets with approximately Gaussian shape moving into 
opposite directions, see also Fig.|5]c. However, there is also 
a region with non-zero Wigner function between these two 
packets that exhibits strong oscillations featuring positive and 
negative values reviling the quantum nature of the coherent 
superposition of the two wave packet s 




5. Conclusions 

In this contribution we have explored ways to visualize quan- 
tum mechanics in phase space by means of the Wigner func- 
tion and the Wigner function flow. The mathematical formal- 
ism of quantum mechanics in phase space in terms of the 
Wigner function resembles classical Hamiltonian mechanics 
of non-interacting classical particles. Thus, teaching quantum 
mechanics in phase space may have some advantages when 
moving from classical to quantum physics. Non-classical fea- 
tures are highlighted in phase space, for example by negative 
values of the phase space distribution or the divergence be- 
tween the direction of the Wigner function flow and tangen- 
tials to the lines of constant Wigner function values of eigen- 
states. 
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FIG. 5: (Color online) Scattering of a wave packed by a potential 
barrier ( 132b in phase space. Potential parameters are Vo = 64 and 
A = 1/8, the wave-packet's initial mean momentum is 8, all values 
in dimensionless units with h = I and m = I. The three upper panels 
show the Wigner function, the position space probability distribution 
and the momentum space probability distribution at different times; 
/ = (part a), t = 1/2 (part b), and f = 1 (part c). Part d illustrates 
the scattering potential i32l . 



